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To reduce noise inside a motor car passenger compartment, two loudspeakers (37 1, 373) are driven by signals de- 
rived from a reference signal x(n) by adaptive filtering carried out by a programmed microprocessor and memory unit (36) 
which adapts the filtering in dependence on error signals e|(n) from four microphones (42,, 422, and 424) distributed in 
the passenger compartment. Reference filtering coefficients are initially determined by analysis of finite impulse responses 
when white noise is acoustically coupled from the loudspeakers (37) to the microphones (42), a white noise generator (48) 
being coupled to the unit (36). The reference signal x(n) is restricted to one or more selected harmonics or subharmonics of 
the fundamental noise frequency by a filter (34) which tracks the selected frequency. The selected frequency may be ob- 
tained from a coil (31) in the ignition circuit of the vehicle. 
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Active vibration control. 



BACKGHOTOn) TO THE INVE^ITION" 

The invention relates to active vibration controls 
As used herein, the term "vibration" includes sound 
or noise, and the invention is particularly concerned with active 
5 noise control • 

In passenger compartments of cars (automobiles), where a 
significant noise component is harmonically related to the rotation 
frequency of a reciprocating engine used to drive the car, the sound 
levels at low frequencies in such enclosures are difficult to 
10 attenuate using conventional passive methods and can give rise to 
subjectively annoying "boom". A method of actively attenuating a 
simple sound field by introducing a single secondary sound source 
driven so that its output is in antiphase with the original ambient 
noise is described in general terms by B. Chaplin in "The Chartered 
15 Mechanical Engineer" of January 1983 1 at pages 41 to 47- Other 

discussions are to be found in an article entitled "Active Attenuation 
of Noise - -The State of the Art" by Glenn E. Warnaka at pages 100 to 
110 in Noise Control Engineering, May-Jiine 1982 and in Internoise 85 
Proceedings, pages 457 to 458 and 46I to 464, and Internoise 84 
20 Proceedings, pages 485 to 488, Particular methods and apparatus are 
also described in British Patent Specification Nos. 1, 577 » 322 and 
2,149,614- 
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SUI-II-IARY OF THE IFmiTIOH 

According to the present invention, an active vibration control 
system for reducing vibration generated by a primary source, is 
•characterised in that at least one reference signal containing 
5 selected harmonics of the said primary soxirce vibration is supplied 
' to means driving a plurality of secondary vibration sources, such 
that vibration energy detected by sensor means operable to sense 
the vibration field established by the primary and secondary sources 
is reduced^ 

10 As used herein, the term "harmonic" includes "subharmonics". 

Preferably the control system is operable in accordance with 
an algorithm which adjusts the outputs from the secondary sources so 
as to substantially reduce a cost function on a time scale comparable 
with the delays associated with the propagation of vibration from the 

15 secondary sources to the sensor means. 

The present invention is particularly concerned with an active 
noise reduction system which can control the sound throughout an 
enclosTore of a car, or at one or a number of "q^uiet zones" within it, 
and which can quickly adapt to changes in the excitation of the sound 

20 field due to changes in, for example, engine load or speed. 
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In order to ensure that the sound produced by the secondary 
vibration sources is of the same frequency as that produced by the 
engine, a signal related to the engine crankshaft rotation rate, for 
example a signal emitted by the engine ignition system, is used to 
generate a reference signal containing a number of sinusoids at harmonics 
(or subharmonics) of the engine crankshaft rotation frequency. 
These are known as engine order frequencies. These sinusoids may 
be obtained using a variety of methods outlined below. Rather than 
attempt to control all harmonics, only a selected set of engine 
order frequencies may, in accordance with the invention, be generated 
as the reference signal. For example, only the firing frequency 
(second engine order in a fovir cylinder car) and its second harmonic 
(fourth engine order) are used if the spectrum of the sound in the 
car is dominated by these components. 

Alternatively, a signal containing all the engine order 
frequencies may be fed to a band pass filter which isolates only the 
particular frequency or frequencies exciting a particular resonance in 
the enclosure interior which could cause a "boom" to be excited. 
The advantage of reducing the number of frequencies fed to the 
filter is that an adaptive filter having fev/er coefficients than 
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would otherwise be the case can be used ThSo m u • •, 

oe usea. This makes implementation 

-or. emcient and allow, fa.ter adaption ti.e. Tha Taata, adaption 
tl^a 13 paxtl=.laPly i.p„„„t automotive applioations, for 

example, in which the active control 

T=ive control system has to adapt surficiently 

quickly to track changes in engine sueed ^H.vv, 

i«xxie speed which may occur on a very 

short timescale. 

Othe. reference aignals be obtained fro. transduoera counted 
on a road wheel hub or the euapenaion ayate» of the car. Such 
reference aisnals „ould contain the harmonica of road wheel 
rotation or road nciae. Transducers placed outside a vehicle ^ 
provide reference ai^s representative of wind noise. Reference 
si^s fro. such sources be broad band (randc). in which case 
these reference si^s ^ be fed directly to adaptive filters as 
described hereinafter. Altamativel., if t.he reference signals are 
periodic (deterministic), control at individual harmonics can be 
exercised as described hereinafter. ' 

'"^ ^P^" i3 a motor vehicle interior, the 

seconder,, sources may be loudspeakers used as the low frequency 
drives of a car audio system. 

Ecamples.of methods in accordance with the invention for 
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generating the reference signals from a signal from an internal 
combustion engine axe now discussed. 
1. Selection of harmonics by filtering, 

A signal is obtained from the primary vibration source which 
5 contains components at all harmonics prevalent in the sound in a 
vehicle powered by an engine. This signal is filtered so as only 
to leave the most important or dominant hermonics. Filtering is 
carried out by a filter whose centre frequency can be controlled 
by an external signal in such a manner that the critical filter 
frequencies have a constant ratio compared with the engine crankshaft 
rotation rate. This can be achieved by using, for example, charge 
coupled devices whose switching frequency is locked to the crankshaft 
rotation frequency, but can also be implemented as a program running 
on a microprocessor, as described hereinafter. 
15 2. Selection of harmonics within a band by fixed filtering. 

A primary source signal rich in harmoxiics is filtered by a 
band pass filter, having a centre frequency fixed at that of a 
pronounced "boom" in the car enclosure and a characteristic such 
that the reference signal only contains the harmonic(s) which are 
particularly exciting the boom. This may be extended such that the 
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filter contains a number of resonances at a number of boom 
frequencies of the car, or even such that the filter has a 
frequency response which models the acoustic response of the 
car interior to the primary excitation. The input signal to the 
5 filter may be a signal from the engine containing all important 
harmonics and may be in the form of a pulse train* 

3. Generation of specific harmonics locked to the engine crankshaft 
rotation frequency. 

This method may be accomplished by using phase lock loops to 

10 generate sinusoidal signals, with frequencies bearing an integer 

relationship to a square wave signal from the engine, which frequencies 
are then atdded together to form the reference signal. Alternatively, 
the signal derived from the engine can be used to control. a number of 
tunable oscillators, each producing a sinusoid at a selected harmoj;y.c 

15 frequency. In one arrangement employing such a bank of tunable 
oscillators, the period of the square wave signal at the engine 
rotation rate is measured with a counter passed to a microprocessor 
which implements a naomber of digital oscillators using difference 
equations of the form 
20 3C^(n) -S(n) +2 cos(rcu )x-|.(n - l) - Xj(n - 2) 

whereto = 2irf-/Nf , I is the order of the harmonic or subharmonic 
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to be generated, is the frequency o f the comter, which counts 
11 pulses during a period, and fg is the sample rate used for the 
difference equation. S (n) is a unit sample sequence to initiate 
the oscillators formed by the difference equations above. The 
5 second and fourth harmonic may be generated for example. 

It is advaintageous if the sample frequency (fg) is derived from 
the counter frequency (f^) by a frequency division circuit for 
example, so that the ratio fs/^c exactly aji integer n;imber. 
An alternative difference equation which may be used to 
10 implement the digital oscillator has the form of a series 
approximation to a trigonometric function, for example 

The variable y is the accumulated phase of the oscillator, given by 

n 

y(n) = Sa^i(n) , 
i=0 

15 ^i(^) ^ where is calculated from the measured period of the 
reference signal, as above, for every new sample (n). The series 
approximation above can be used for y(n) in the range -'^/2^ y(nX Tr/2. 
For values of y(n) outside this range the symmetry properties of the 
cosine waveform are utilised, until y(n) >1T. V/hen y(n) , the 

20 natural overflow properties of the two's complement number 
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representation is used to allow y(n) to "'.vrrap-around" to 
y(^) "TT and the series expansions and symmetry properties 
discussed above are again used. In this way the digital representation 
of y(n). is kept within the range - '^-^ y(n)-<'Tr, and Xj(n) is within 
5 the range - 1. 

The calculations of the coefficients used in the difference 
equations forniing the digital oscillators, together with the 
difference equations themselves may be implemented on a dedicated 
processor, or may form part of the program which also implements 
10 a controller that generates the outputs used to drive the secondary 
sources from the reference signals described above. 

The controller is designed to be adaptive so as to quickly 
track changes in engine speed and load. The outputs of the secondary 
sources are adaptively controlled so that some measurable cost 

15 function is minimised. This cost function would typically be the 
sum of the mean square outputs from a number of microphones in the 
enclosed space. The controller can be implemented as a digital 
adaptive FIR filter, using the basic update algorithm described by 
S. Elliott and P. llelson in "Electronics Letters", at pp. 979-981, 

20 1985. A nxamber of additions must be made to this basic algorithm, 
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however, to enable it to work quickly and emciently in this 
particular example. THe basic algorithm, referred to hereinafter 
as the stochastic gradient algorithn, is presented below, in order 
to highlight the necessary alterations. If the i-th coefficient 
y . Of the adaptive filter driving the .'th secondary source at sample 
number n is w^^Cn). then each of these coefficients should be 
adjusted at every sample according to the equation 

w^i(n + 1) = w^i(n) +o<Str (n . i) 

where oc 13 a convergence coefficient, e (n) fH^ . 

x^ieiiG, e^^n; is the sampled output 

from the Z -th sensor and r Jn) is a sequence formed by filtering 
the reference signal discussed above ( x(n), say) with a digital- 
filter which models the response of the^ -th sensor to excitation 
Of the m'th secondary source. These digital filters, generating 
-ach r^ Jn). have only two coefficients in the implementation 
described in the "Electronics Letters'- article, since control at 
only a single, fixed, frequency was being attempted. In the present 
example, however, the digital filters ^st model the relevant response 
over a range of frequencies, governed by the frequency range which 
the active system is attempting to control. It has been found that 
under certain circumstances a digital filter need only model the 
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overall delay in the response to ensure the stability of 'an adaptive 
filter. It is more common, however, to have the digital filters 
incorporate a delay and then some reverberant response. This may 
be implemented using either digital FIR or IIR filters whose 
5 coefficients are adjusted adaptively during an initialisation phase, 
so as to accurately match the desired responses. It is also possible 
to continue this initial adaption process during the operation of the 
•active control system by feeding training signals to each secondary 
source which are suitably uncorrelated with each other and with the 
I primary excitation. This may be necessary to track changes in the 
acoustic response of the enclosure. Alternatively, if the chaage is 
due to some well-defined cause, such as a passenger sitting down or 
a car window being opened, this change may be detected with 
mechanical transducers and the information used to switch between a 
cariety of filters modelling the response of the enclosure under a 
variety of conditions. 

Another important consideration concerning the use of adaptive 
algorithm in this example concerns the effect of unwanted, low level, 
harmonic or subharmonic components in the reference signal. Suppose' 
that the method used to generate the reference signal, as described 
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above, is designed to produce only I frequencies. 3ven though 
there may only be I desired harndnics in such a system, there will 
in practice also be a number of other harmonics or subharmonic 
frequencies at low level, because of the finite cut-off rate of 
? the filters, for example. These components aay also be being 
generated by the primary source and therefore present in the 
enclosure aiod hence at the outputs of sensors such as microphones 
and thus the adaptive algorithm will attempt to cancel them by 
enormously amplifying, the low level, spurious harmonic reference 
signals. This can cause numerical overflow problems in the adaptive 
filter coefficients. This may be prevented in a num^ber of ways: 

(1) The use of only 2.1 coefficients in the adaptive filter. 

(2) The deliberate injection of random noise into the 
. reference signal. 

15 . (3) The insertion of a "leak" into the algorithm so that, in 

the equation above, the past coefficient value, w . (n) 
is multiplied by a factor close to but not equal to unity 
before being updated. 
(4) The -addition of an extra term in the update equation which 
minimises a cost function involving "effort" as well as 
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"error", as described in ISVR Technical Report llo. 136, 

1985 1 published by the Institute of Sound and Vibration 

Reseaxch, University of Southanipton. 

A n\imber of other suiaptive algorithms may also be implemented 

5 to adjust the coefficients of the digital filters in the controller 

driving the secondary sources. These alternative algorithms are best 

described in matrix form. 

Assuming the availability of a sampled reference signal, .x(n), 

v/hich is correlated with the output of the primary source, but is 

10 ixnaffected by the action of the secondary sources. The output to 

the m'th secondairy soiirces, yjjj(n), may be obtained by passing this 

reference signal through a digital filter v/hose i*th coefficient is 

Wjj^(n) at the n*th sample, so that 

I-l 

yia(n) - w^{n)y:{n - i) 

i=0' • ' . 

15 The sampled output from theX 'th error sensor, e^(n), is equal 

to the sum of the contributions from the primarrr source, d£ (n), and 

each of the secondary sources. The response of the path between the 

m'th secondary source and^ *th error sensor is modelled as a J'th 



20 



order PIR* filter with coefficients c^rQ-j so that 

M J-1 ' ^ 
e£ (n) = di(n) +2 S c^ajymC^ " j) 
m=l j=0 
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Therefore 

m=l jTo '-^"J ^n^i^^-^^^^^-i-j) 
in order to obtain a ^trix expression Tor the error surface 
we .ust no. .a.e the assumption that the filter coefficients in 

' - ' -e controXXer onX. adapts 

ver. sXowX. compared to the time scaXe of the response of the 
syste. to he controXXed. ^Then w,,(n-,) ^ 

e^(n) = d^(n) + 2 S ^ w ^ • ^ 

If v/e Xet 

J-X 

''^r.M = _S c^^.x(n-o) 

Where the sequences r,Jn) for each ^ and . a.e caXXed the fiXtered • 
reference signaXs, then 

Dl=l 1=0 

where 

^ ^jCn) ... ^^^^^^^ r£„Cn-I+X)] 

- =" "XO WnT , -I 

"30 • • • 

So if 

1%) =[di(n) d2(n) ... d^Cn)] 



MI-X 
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10 



then 

e,(n) = d(n) + R(n)w 

where 

R'^(n) =[r3^ (n) (n) ... (n)] 

If the cost function is written as 
L 

J = E( S ei^in) ) = Ei:/(n) e(n)J 
-C =1 

where B is the expectation operator, then 

J = B(dT(n)d(n) ) + 2wTE(RT(n)d(n) ) + w^E(RT(n)R(n))w 
Using the standard theory of matrix quadratic foras, the 
minimum of J, Jmin, is obtained with 

iiopt = - -E(R'^(n)R(n))-lE(RT(n)d(n)) 

giving 

; Jfflin = c - i^A'-^b 

= E(dT(n)d(n)) - SCd^(n)R(n))E(R-(n)R(n))-lE(RT(n)d(n)) 
^5 The true gradient may be written 

}i " 2E(RT(n)d(n) + ^(n) R(n) w(n)) 
aiJ = 2E(RT(n) e(n)) 

The time domain steepest descent algorithm may thus be written as 
'■^^1"^- 2^E(R'I^(n) e(n)) 
20 In a practical implementation the true expectation could be 
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approximated by an WA or AR averaging process. Alternatively, the 
. instantaneous gradient could be used to update each filter 

coefficient every sample, as in the time domain "stochastic gradient" 
algorithm 

i w(n+l) = w(n) - ^^(ti) e(n), 

which is the matrix representation of the algorithm described above. 
It is clear from this fomulation that the stability and convergence 
properties of this algorithm, in the limit of slow adaption, are 
governed by the eigenvalue spread of the matrix E{R^{n) R(n)). This, 
matrix depends only on the response of the system to be controlled, 
the positioning of the sources and sensors vnLthin that system and 
the spectral properties of the reference signal, x(n). it is 
possible that an unfortunate placing of these sources and sensors 
could-cause this matrix to become-ill conditioned, so that it has a 
large eigenvalue spread. This would create very slow "modes" in the 
convergence properties of such an algorithm. 

This problem could be removed by using a ITewton's method 
algorithm, the exact form of which may be written as 
ifk+1 = Wfc - 2^E(rRT(n) R(n)J rh(R^(n) e(n)) 
Again, various types of averaging could be used to give a 
practical approximation' to the expectation operator. It should be 
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noted however that the tirie-independent matrix E(R^(n) R(n)) 
depends only on the response of the system to be controlled amd 
on the reference signal, and these are assuaed to be knovm and 
stationary. This suggests a variety of stochastic Newton's method 
5 ("SIM") algorithms* The most obvious of these uses a modified or 
"normalised" set of reference signals (n) , such that 

The computation of each of these reference signals will take 
somewhat longer than for R(n) alone since none of the elements of 
10 Q (n) are necessarily time delayed versions of any other elements. 
The complete SNM algorithm, again using instantaneous versions of 
E(R'^(n) je(n)), becomes 

BRIEF DSSCRIPTIQN OF THE DRA^/HINGS 
15 anbodiments of the invention will now be described, by way of 

example only, with reference to the accompanying drawings, wherein 

Figure 1 is a block schematic diagram of an active noise control 
system associated with an enclosed space, 

Figures 2(a'), (b) and (c) are graphical representations of the 
20 behaviour of an element of the system, 
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Figure 3 is a block diagram of one fo^ of reference signal 
generator, 

Pi^e 4 U a .locic diagram of another for» of reference ei^x' 

generator, 

5 n^re 5 is a bloc, diasra. of a circ.it which inoorporatee a 

microprocessor, 

^ n.'ure 6 is a bloC diasra. illustrating hcv, t„o reference 
signals are combined, 

Figure 7 is a flow chart illustratinc. r,^^^ 

iixustrating processing in a particular 

10 embodiment of the invention, . 

a^e 8 is a soheoatic dia«ra« illustrating a heterodyne and 
averaging .ethod of obtaining inphase and quadrature exponents of 
an error sequence in another embodiMent of the invention, ' 

Figure 9 illustrates application of the invention to non-noise 
15 vibration control, and. 

Figure 10 illustrates a modification 

moaiixcation of the arrangement shown 

in Plgure 9, 



20 
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pgrAILSD DESCRIPTION OF THE PR2FI51RED SIZ?-0DES::TS 

In Figure 1 an enclosure 10, which is the interior of the 
passenger or driver compartment of an internal combustion engine 
driven vehicle, in this example, a motor car 100, is represented 
5 schematically together with an active sound control system 1 

according to the invention. In this example, the system 1 employs 
two secondary sound sources 11, comprising two low frequency 
- loudspeakers of a stereo audio system fitted to the car, arid 
three acoustic sensors, comprising microphones 12, The loud- 
10 speakers 11 are driven by a controller circuit 13 which comprises 
a pair of adaptive filters 14. Each adaptive filter 14 drives a 
respective one of the loudspeakers 11 with an output signal 3 which 
the filter 14 produces as a-result of its 'action on a reference 
signal 4 supplied thereto by a reference signal generator 15* The 
15 reference signal 4 is generated by the generator 15 from an input 
signal 16 which is periodic at the crankshaft rotation 
rate of the internal combustion engine 2. 

The signal generator 15 may comprise a tracking filter. 
The purpose of the outputs from the loudspeakers 11 driven by 
20 the controller 13 is to reduce the sound vibration field established 

by the primary and secondary sources, experienced within the enclosxire 10 • 
Since the primary source (engine 2) of the noise to be reduced is 
periodic, the reference signal 4 generated by the generator 15 is, 
in accordance with the invention, arranged to' contain one or more 
25 sinusoidal components at harmonics (or subharmonics) of the crankshaft 
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routioa rate of the engine 2. «,e adaptive filter. 14 axe adjusted 
automatically output .ignala 5 fro. the sensor-eiorophones 12. 
oorresponding adjustment being .ade simultaneously to the outputs 
of the loudspeaj^ers 11 so as to substantially miniMsa a cost function 
. on a ti.e scale comparable with the delays associated with the 
propagation of sound .vibrations fro. the loudspealcers 11 to the 
-iorophones 12. The cos. i^^ction .ay comprise the sul of the mean 
square outputs of the microphones 12. 

A decision is made beforehand as to wha+ >, >^ • 

-licuia as ^o what harmonics are to be 

selected; a decision which may vary from car to cax. 

The control system 1 does not employ a stored solution. Instead, 
xt makes use of a plurality of closed loops, each loop comprising a 
Microphone 12, the controller 13. and a loudspeaicer.il, whereby 
Signals from the microphone 12 are used to adapt the filters I4 
controlling the loudspeaker 11. whi.h has an influence on the output 
Of the microphone as a result of acoustic response within ^he 
enclosure 10. 

&=h loop is accounted for by part of an algorithm which adjusts 
the outputs Of the loudspeakers U as aforesaid, the algorithm being 
Of the form 

w(n^.l) = w(n) R'^(n) e(n). 

It will be noted that the loudspeakers 11 and the microphones 12 
are distributed in the enclosure 10 in spaced relationship.. The 
distribution, which varies from car to car. is adjusted to get 
substantial sound reductions throughout the enclosure 10. ' 

It will also be noted that .the system 1 employs as . many closed 
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loops as the niimber of sensors (12) multiplied by the number of 
secondary sources (ll). 

The manner in which each secondary source 11 effects every one 
of the sensors 12 is reflected in the algorithms referred to herein. 
5 Furthermore, the system 1 employs more sensors (l2) than secondary 

sources (ll), vhereby a controlled reduction of primary source vibration 
is achieved. This contrasts with presently- known systems employing 
"the same n\mber of secondary sources as sensors, whereby near perfect 
cancellation can be achieved at the sensor locations but vibration 

10 levels away from these locations are increased. 

Since the reference signal 4 contains harmonics of the input signal 
16, which signal is periodic at the engine crankshaft rotation rate, the 
reference signal 4 contains engine order frequencies. The signal 
generator 15 is arranged to select engine order frequencies that ensure 

15 that the sovmd produced by the loudspeakers 11 is of the same frequency 
or frequencies as the sound produced in the enclosure 10 by the engine 2, 
even during changes in engine conditions such as load or speed. The 
number of engine order frequencies in the reference signal 4 is 
restricted so that the adaptive filters 14 have a relatively small 

20 number of coefficients and can therefore adapt quickly. 

All the coefficients- are constantly being updated by the system, on 
a sample by sainple basis. Thus there is no waiting for a final response 
before making another adjustment. The sample time is only a small 
fraction of the fundamental frequency of the primary source 2. 

25 The input' signal 16 can be obtained from another moving part of the 

engine or part of the ignition circuitry; for example, 

Fig-ure 2(a) illustrates graphically the response of the reference 
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generator 15 when in th. for. of a Wicins ruter. l^at 
is to .a., a filter having a centre f.e,uenoy .o oontrolled that 
the filter output frequencies have a constant ratio to the .cinant 

input frequencies, so that in Figure 2(bl) +ho r 

xgure ^^a;, the frequency fo is N x 

' (engine crankshaft rotation rateV an^ +k 

^ the reference signal 4 

(Slaure l) contains only the fiT.«?+ w h 

y ^ne xirst N harnonxcs of the engine 

- a voltage pulse train as represented hy 3, 
xs the periodic ti.e of the engine rotation rate, the first eight 

haraonics of the engine rotation rate ar. r... ^ ■ 

-Lon rate are present m the reference 

aisnal 4. The spectrum, by Fourier analysis of ^h. . r 

°^ reference signal 

4 xs then illustrated by Figure 9/'^ ^ • , • , 

7 ^^Sure 2(c), m whxch A is amplitude. With 

the response of Figure 2(a), only the fiT.<=+ • u 

\ J, onxy tile first sxx harmonics would be 

usable. 

The trao:.an. filter comprising the signal generator, can he in ' 
the for. Of ch^ge coupled devices having a switching frequency 
locked to the engine crankshaft rotation rate. 
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An alternative for=i of reference signal generator, which 
employs a plurality of tracking band pass filters is illustrated 
by Figure 3 in which the input si^l l6, a square wave at, for 
example, 128 ti.es the engine drive shaft rotation rate, is divided 
> first V 32 and then ,y 2. division by 32 is achieved by a divider 6 
Which produces a s.uare wave si^l j four ti.es the engine cranio 
shaft rotation rate, which signal is supplied to a bandpass filter 17 
The filter 17 has a centre freouancy f^ which is arranged to 

track the fundamental freauencv r.r 

-Lxequency or the square wave signal 7 

supplied thereto. The further division ^ 2 is achieved e divider 

8 Which prcducee a s^uaxe wave eignal 3 at twice the engine drive 

Shaft rotation rate. The ei«^l 5 is then supplied to a handpass 

filter 18 having a centre frequency and which is axra-nged to track 

the fundamental frequency of the square wave signal 9 supplied 

thereto. The bandpass filter<? 17 or,^ lej 

F Hirers 17 and 18 produce respectively 

Sinusoidal cut^t signals 7a, 5a at f, and f, which are linearly 
sunned at an adder 1, to produce the required reference signal 4. 

Itether dividers and tracking handpass filters can, of course, he 
incorporated in the circuit of n^ure , so that the reference signal 4 
contains the desired set of engine order frecuencies. 

Another firm of reference signal generator may comprise a fixed 
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frequency filtering circuit which selects harLionics and/or 
subhamonics from axi input signal rich in the harmonics of the 
engine crankshaft rotation rate or firing rate. The filtering 
circuit may comprise a bandpass filter having a centre frequency* 
5 fixed at the frequency of a pronoxinced resonance excited in the 
enclosure 10 (Figure l) by the engine 2 or other primary source 
of vibration. For example, a bandpass filter may be arranged to 
have a response that models the acoustic response of a motor car 
passenger compartment to the engine. 
0 A further form of reference signal generator may comprise a 

plurality of phase lock loops used to generate sinusoidal signals 
having respective frequencies with integer relationships- to a 
square v/ave input signal from the engine 2 or other prijnary vibration 
source. The sinusoidal signals may then be added together to form 
5 the required reference signal. Thus a reference signal 4 comprising 
specific harmonics and/or harmonics locked to the primary source 
fundamental, such as engine crankshaft rotation rate, is generated. 
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An alternative generator for such a reference signal ig 
illustrated in Figure 4 in which a square wave 20 at the primary 
source fundamental is used to control a plurality of tunable digital 
oscillators 25, 26 each producing a sinusoidal signal at a chosen 
i harmonic oi; subharmonic frequency to be added at the adder 19 which 
produce the reference signal 4 by simple addition of the sinusoidal 
signals. 

In the signal generator 27 of Figure 4. the square wave signal 
20, which is at an engine crankshaft rotation rate, is supplied to 
a bistable circuit 21 which divides the rate signal by two and thereby 
produces a pulse train signal 20a in which the duration of each pulse 
is equal to the prevailing periodic time of the square wave signal 20. 
This periodic time is then measured by a counter 22 which is enabled 
throughout the duration -.of each positive pulse from, the bistable 
circuit 21 and "counts clock pulses supplied by a clock pulse 
generator 23. The clock pulses axe generated at a fixed, suitably 
high rate fc_. 

The contents of the. counter 22 are read at the end of each 
positive pulse from the bistable circuit 21 by a trigonometric function 
generator 24. The generator 24 is triggered by the trailing pulse of 
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each positive pulse from the circuit 21 and generates two digital 
outputs representing respectively cos (2 q) and cos (4 q) where q 
is given by 

" 0 = 2lTf^/Nf^ 

5 in which M is the number of clock pulses counted by the counter 22 
in the duration of one positive pulse from the circuit 21, and f„ 
is a sample rate used in the two digital oscillators 25 and 26 which 
receive respectively the digital outputs cos(2<^) and cos (4^) from 
the function generator 24- The digital sinusoidal outputs from the 

10 tv/o oscillators 25 and 26 are superposed by a digital adder 19"^ 
which supplies the reference signal 4 as a digital signal. The 
trigonometric fiinction generator 24, oscillators 25 and 26 and adder 
19^ can be implemented by a microprocessor with a suitable program • 
The microprocessor can be used to produce the reference signal 4 in 

15 the form :- 

XjCn) = S (n) + 2 cos(l'^)xj(n - l) - Xj(n - 2) 
where I is the order of the harmonic or subharmonic being generated, 
S (n) is a unit sample sequence which initiates the simulation of 
the oscillator, and n is the sample number. 

20 
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5 Figure 5 represents in block form an active sound control 

syS^tem 30 for reducing the level of engine generated noise in the 
passenger compartment of a motor car. The car is provided with an 
ignition circuit including a low tension coil 31 from which a 
voltage signal 32 at the firing rate of the engine is taken and 

10 supplied to a waveform shaper 33 which in response thereto produces 
a pulse train at the engine firing rate. It is assxamed in the 
present example that the engine firing rate is twice the engine 
crankshaft rotation rate fQ. Thus the shaper 33 provides a signal 
having a fundamental frequency which is a single harmonic, (2fQ), of 

15 the crankshaft rate. A reference signal generator is provided 

in the form of a proprietory tracking filter 34 j manufactured by 
Bruel and Kjaer under type nxmber 1623. The tracking filter 34 
receives the output of the shaper 33 as an input signal and as a 
trigger signal and produces a sinusoidal output signal at the selected 

20 harmonic 2f^. This sinusoidal signal is sampled with an analog to 
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digital converter 35 to produce a reference sequence x(n) of 
digitised ssunples which are supplied as data to a processor and 
memory unit 36. 

Mounted within the motor car passenger compartment 10 (not 
5 shown in Figure 5) are two loudspeakers 37 i and 372 1 ^^^^^ ^ 

positions normally used for car stereo reproduction.. The loudspeaker; 
37-^, 572 are driven by a multiplexer 38 through respective low pass 
filters 39 and output amplifiers 40. The filters 39 have a cut 
off frequency of 46O Hz and are provided to prevent aliasing. The 

10 multiplexer 38, which contains sample and hold circuits for each 

output, is controlled by the processor and memory unit 36, through 
a control line 55, and receives a single input signal 57 from a 
digital to analog converter 4I. The purpose of the loudspeakers 
37^ and 372 is to generate, in the passenger compartment, audio 

15 waves that will cancel those set up directly by mechanical 

transmission from the engine to the compartment. The digital to 
analog converter 4I is supplied by the processor and memory ^xn±t 36 
with output data 58 which consists of two interleaved sequences of 
digitised samples y^(n) and y^Cn). The data 58 is converted by the 

20 converter 4I into interleaved sequences of analog samples and 
separated into respective sequences by the multiplexer 38 for 
application to the low pass 
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10 



filters 39. Thus, in effect, the loudspeaker 57i is driven by 
the sequence y-j^(n) and the loudspeaker driven by the sequence 

y^Cn). In Figure 5 each sequence of data 58 is represented by the 
expression y^Cn), so that in this example m may be 1 or 2. 

In order to ensure that the acoustic outputs from the loudspeakers 
37]^ and 372 correct phase and amplitude to effect 

cancellation of the engine noise, error signals are picked up from 
the passenger compartment and utilised by the processor and memory 
unit 36. Acoustic error signals, if present, are sensed by four 
micrpphones A2-^i 422 > 42^ and 42^, which axe placed respectively 
either side of a driver headrest and a passenger headrest, there 
being only two seats in the compartment in the present example. The 
electrical outputs from the microphones 42^ etc, are respectively 
amplified by amplifiers. 45 and passed through low pass filters 44 
.5 to a four-input multiplexer 45 which supplies a single analog output 
to an analog to digital converter 46. The filters 44 are provided 
to prevent aliasing and have a cut-off frequency of 46O Hz. 

The multiplexer 45 is controlled by the processor unit 36 
by way of control line 56. 

)0 The multiplexer 45 and the converter 46 convert the four 

filtered microphone outputs into a data stream 59 comprising four 
interleaved sequences of digitised samples e^Cn), SgCn), e^(n) and 
e4(n), which correspond respectively to the filtered outputs of the 
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.icrophcnea ,Z,. 42„ 42, and 42,. In K^.e 5. each sequence is 
represented ^ . („) „ ^^^^^^ ^ ^_ ^ ^ 

Ihe processor and .e.or. unit 36 receives a s,.axe wave si«naifco 
at 1.2 kilohertz from a samnlo 

sample rate oscillator 47 which detex^inea 

the rate at which the converters 33, 41 ^ 46 convert samples aM 
the rrame duration or processing carried out the unit 36. .hus 
xn the present example, the unit 36 completes its processing frame 
Within 833 billiseconds. , crystal clock oscillator 6I with a 
frequency of 10 Megahertz is included in the unit 36 



The unit 36 simulates two adaptive filt^r-^ v, u • • 

udpgive I liters, each having two 

coefficients, so that 

ym(n) = w^Qx(n) + w^^x(n - i) 
describes the relation between the output sequence .,(n) to a 
loudspeaker and the reference signal .(„), „,e.e the coemcients 
are and w^^. Hence with the two loudspeakers 37i and 37^ ,. 
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y^Cn) = W2Qx(n) + W2ix{n - l). 

The values of the coefficients w _ and w . are calculated by 

mO ml 

the unit 36 from the relationship 

4 

^mi^^ + 1) = ^mi^""^ "^^^^ ®-d^^)^£ni^^ " 

5 

in which oc is a fixed convergence coefficient, r^^(n - i) is a value 
of a filtered reference signal r- , and i = 0 or 1. 

The filtered reference signal r^^ is a seq,uence formed by 
filtering the reference signal x(n) with a filter that models the 

10 . effect of the acoustic coupling between the m"^^ loudspeaker and the 

^ microphone. The unit 36 simulates this filtering as digital PIR 
(Finite Impulse Response) filtering. Coefficients for the digital PIR 
filtering are adjusted adaptively during an initialisation program in 
which a white noise generator 48 is energised. 

15 In the initialisation program, a white noise signal is generated 

by the generator 48, which is then filtered by a low pass filter 49 
to prevent aliasing, the filter 49 having a cut-off frequency of 
460 Hz. The signal is .subsequently sampled and converted by an einalog 
to digital converter 50. The digital output of the converter 50 is 

20 used to drive the loudspeakers 31 1 and 37 2^ °y ^'^^.y of the processor 
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and memory unit 56, and the resultinc digital input to the unit 36 

from the microphones A2-^, i>,2^, 42^ and 42^ is used to determine 

the values of reference filter coefficients 0^^^^^ where j => 0,....,34. 

The unit 36 performs a 35 coefficient PIR modelling of the impulse 

5 response between the m*^ loudspeaker and the -t*^ microphone at the 

j sample. Such modelling is described in "Adaptive Signal 

Processing" by 3..V/idrow and S.D. Stearns, published in 1985 by 

Prentice Hall. 

The filtered reference sequence is then given by :- 

34 

j=0 

The operation of the unit 36 is such that, having obtained the error 

samples e^(n) and the filtered reference si^al r- (n), each adantive 

cm 

filter coefficient w^^^ for each output y^(n) is updated by a quantity 
proportional to the sum of the 'computed products of e£(n) and 

15 ^^in(^ " ^) accordance with the equation : - 

\ 4 ! i 

^mi(^ + 1) = *^mi(^) +^ 2 e.(n)- ^Cn - i) 

^=1 

The new set of adaptive coefficients Wj^^ is then stored and used to 
filter the next sample of .the reference signal, x(n +1). 

The unit 36 includes RAM for temporary storage and computation, 
20 and EPROM for program storage. Calculated coefficients w^^j^ and c^^^j , 
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. and reference sequences r^,(n) are held in HAir. The convergence 
coerricient is entered at a set of .anuall, operable switches (not 
shown) . 

preferably the unit 36 includes a Texas Instruments TMS 32010 
Microprocessor. The input signal rate fro. the ignition circuit, 

waveform shaper 53 is a monostable circui- f . . v 

aoxe cxrcuio triggered by the leading 

edge of the input signal to produce duIs^, 

proQuce pulses of a constant width which 

1= reutive to t..e .a.pZe perio. .et t.e sample «te „r z a 

kilohert3. The low p... r,!,,,^ 44 ^ 49 are active mter 
modules supplied by Kemo United uaier No. 1431/1. 

only a a^l amount or separate additional is required with 
the ahove-entioned «s „oio microprocessor, which has cohsiderahle 
internal H^: an, operates as described in the T,E ,2010 .sers. Guide 
published in 1,33 by Te:=as Instrueents Inc. Bata buses between the 
unit 3« and the converters 35, 41, 46 and 50 are 12 bit buses, other 
^uses and lines required Tor synchronisation and control are omitted 
ror Clarity. It will be noted tbat the values o. the coemcients 
"mi ^ c^mj can be initially set to zero. 



10 



15 



!0 



WO 88/02912 

PCT/GB87/00706 

- 33 - 

A family of algorithms which are alternatives to the Smi 
algorithm uses a bank of adaptive digital filters working in 
parallel for each secondary source. Ekch individual filter is fed 
by a reference signal containing a subset of the harmonics or 
i subharmonics to be controlled. For example, Figure 6 of the 

accompanying drawings shows two parallel FIR filters 70 each fed by 
pure tone reference signals 71, at the second aiad fourth engine 
order frequencies in this case. The outputs 72 of these filters 
are added together by an adder 73 to form an output 74 to the 
secondary source. Each of the parallel filters 70 may be updated 
by any of the. algorithms discussed above; For example, the said 
stochastic gradient algorithm may be modified' so that 

WjmiCn * 1) = wi^(n) | e^(.^) ^^^^^^ ^ 
where w^ is the i'th coefficient of thellR filter fed from the 
I'th harmonic of the engine, driving the m'th secondary source.. 

The advantage of such an algorithm is that each harmonic ' 
frequency is controlled independently and the convergence of an 
harmonic does not couple with the convergence of any other harmonic, 
as is the case when a 21 coefficient filter is used to filter I 
harmonics simultaneously. . In this case I filters each with 2 
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coefficients could be used to fin-^-^ t u 

used to filter I harnorucs individually, and 

their response combined afterward=! tv,^ ^ • j o. 

dj.-erwaras. The disadvantage of this 

al.o.ith„ that a nite... „rere„=e si^.l „.^, 

for each source (.), sensor M aM harmonic (I), to each 

. mother approach to controUins a number or harmonics is to ta^e 
the courier transfer, of each of the error 3i.-nals an. to update a 
set of coefficients controlling each harmonic of each secondary 
output inaependentl,. The outputs of each of these filters are then 

combined together, for each secondary sou-ce via =r, ■ 

source, via an inverse Fourier 

transfer,, to generate the output wavefor. for this source, as 

indicated in Figure 7. 

- For a single harmonic in the fre.uenc. domain, the complex value 
of the -th error signal will be given by 
E * A 2 

''here A, is the value If C^th^n: active control, '„ is the 
complex amplitude of the voltage to the .Hh secondary source aM 
is the coiple, transfer function between the^-th sensor and 

m'th source at the frequencv of th^ v, • 

^equency of the harmonic of interest. 

Ti.±s may be expressed in matrix form as :- 
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i; - (El S2 ... 



^"^^ (^1 ^2 ... v/,p 

^ = ['^11 h2 •••• C^m) 



( 



21 ) 



( : 



) 



The cost function in this case .ay be written as J - where the 
superscript H denotes the complex conjugate of the transposed vector 
or matrix. Therefore 

J = A^A + W^c^A + A^c W -K W^C^^ 

Note that 

so that the steepest descent algorithm may be -bitten 

Where Wjc and ^ are the complex values of the filter response and 
error output respectively at the k'th iteration. 

This algorithzn does not seem to appeax in the active control 
literature. This is probably because the .ewton-s method algorithm. 
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below, is no more difficult to implement after initialisation to 
form the matrix needed to prenultiply Hov/ever, in some cases 

it may be that the matrix £ changes with time and a separate 
"identification" algorithm is used in parallel ^/d-th the adaptive 
5 control algorithm to track these changes. In such cases the 

steepest descent algorithm may be considerably more computationally 
efficient to implement than that below. 

The frequency demain version of the ITewton's method algorithm 
may be 'bitten 

A special case of this algorithm does appear in the literature, in 
which the number of error sensors is equal to the nximber of secondary 
sovirces (L = M) , so. that £ is a square matrix, and the algorithm 
reduces to 

This algorithm has been presented by Pierce (1985, David W. Taylor 
Maval Ship Research and Development Center Report No. 85/047. An 
algorithm for active adaptive control of periodic interface). If the 
convergence coefficient,/^, is set equal to one half, the algorithm 
3 also reduces to the iterative matrix algorithm described in 

White and Cooper (1984, Applied Acoustics 17, 99-109. "An adaptive 
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controxu. fo. .uZU.„i...e activ. control..). Sea .I.o U.K. 
Patent Specification Ho. 2,122,052A. 

in practice the .ost advantageous algorithu „iU probacy be 

derived from a judicious mixture of tl»= <■ 

i=ure or time and frequency domain concepts. 

I'or exampie, one implementation of a Pourier transformer operating on 

eao. error sequence e^ (n) to produce the inphase and quadrature 

frequency components of e^ (n) at I„ is illustrated in lU.-ure 6 

Wherein integrators SO and multipliers 81 are used. The slowly 

varying outputs of this circuit represent th. r„i . • 

'l"'-^^"^ the real and imaginary oarts 

of the frequency domain signal defined above, so these signal^ 
can be used with any of the frequency domain algorithm, discussed 
above to update a ban. of adaptive filters driven hy this frequency 
component drivir.g each secondary source, as in Figure 6. 

il^e 9 illustrates application of the invention to non-noise, 
■5 I.e. neohanioal, vibration control. 

The example illustrated by Pi^^e 9 comprises a modification of 
the Blg^e 1 arrangement, wherein the microphone sensors have been 
repucsd by accelerometers 90 and the loudspeakers by mechanical 
Vibrators 91. The accelerometers 90 and vibrators ,1 are mounted 
5 on surface portions of the enclosure 10. 
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In another modification, illustrated by Figure 10, a sensing 
combination of microphones 12 and accelerometers 90 is used, and/or 
a source combination of loudspeakers 11 and vibrators 91 • 

A combination of Figure 10 need not be coincident; the 
5 components thereof could instead be spaced from each other. 
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CLAIMS 

1. An active vibration control system for reducing vibration 
generated by a primary source is characterised in that at least 
one reference signal (4) containing at least one selected harmonics 
o:^ the said primary source (2) of vibration is supplied to means 

5 driving a plurality of secondary vibration soxirces (ll), such that 
the vibration energy detected by sensor means (12) operable to 
sense the vibration field established by the primary and secondary 
so\irces is reduced. 

2. A system according to Claim 1, characterised in that the 
10 reference signal (4) is generated by filtering a periodic input 

..signal (16) having its fundamental frequency (9) locked to a 
predomin^jit frequency (fg) of "the exterior soijrce. 

3« A system according to Claim 2, characterised in that filtering 
is effected by at least one tracking filter (17» IS). 
15 4» A system according to Claim 2, characterised in that the 

reference signal (4) is generated by at least one tunable oscillator 
(25) whose frequency is controlled by a signal (20), the primary 
source fundamental. 

5» A system according to any preceding claim, characterised in that 
20 the means (56) driving the secondary sources (57> 522) comprise a 
digital processor with data and program memory. 

6. A system as claimed in any preceding claim.- characterised in that 
means (70) are provided whereby a plurality of reference signals (71) > 
each containing a single harmonic, are filtered by a bank of filters (70) 
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which are adjustable independently eind whose outputs (72) are 
coabined (74) to form an output to the secondary sources. 
7. A system according to any preceding claim, characterised in that 
it is adapted to operate in accordance with an algorithm which adjusts 
5 the outputs from the secondary sources so as to substantially minimise 
a cost function on a time scale comparable with the delays associated 
with the propagation of vibration from the secondary soiirces (ll) to 
the sensor means (12), 



8, A system as claimed in Claim 7, 


characterised 


in 


that 


the 


algorithm is of the form 










w(n+l) = w(n) -p^^(n) 
9* A system as claimed in Claim 7, 


e(n). 

characterised 


in 


that 


the 


algorithm is of the form 










iin:+l = ^ - 2^§^(n) e(: 
10* A system as claimed in Claim 7, 


a). 

c har ac t er i s ed 


in 


that 


the 


algorithm is of the form :- 










W = W - 2>uC^E 
11. A system as claimed in Claim 7, 


c har ac t er i s ed 


in 


that 


the 



algorithm is of the form :- 

12. A system according to any preceding claim, characterised in that 
it makes use of a plurality of closed loops (12, 13, ll), each 
comprising one of the secondary vibration sources (ll), signal 
adaptive control means (14) and one of the sensors (12). 
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13. A system according to any preceding claim, characterised in 
that the vibration controlled comprises soiand or noise. 
14- A motor car provided with an active vibration control system 
as claimed in any one or Claims 1 to 13 . 
5 IS* A system for reducing vibration, substantially as hereinbefore 
described with reference to the accompanying drawings. 



wo 88/02912 



PCT/GB87/00706 



1 17 

FIG.1 



100 



5 



1 



10 




ENGINE 



15 



/■2 /13 



16 

JU 



REFERENCE 

SIGNAL 
GENERATOR 



ADAPTIVE 
FILTER 



ADAPTIVE 
FILTER 



3- 



3 



1 



SUBSTITUTE SHEET 



J 



wo 88/02912 



PCT/GB87/00706 



2/7 



FIG. 2 




fl — JL^t (b) 



8T 9T 16T17T 



(c) 



123i»5 678 



Nf 



'/fc 



r 



15 



FIG. 3 



t32 



6 



^17 



-TLTLrLn 



7 



7 




7 



8- 



19 



r2 



9 



7 




•18 




9a 



SUBSTITUTE SHEET 



wo 88/02912 



PCT/GB87/00706 



3/7 



LTL- 
i 



20 



DIVIDER 
-2 



20a 



21 



FIG.Z. 



COUNTER 



TRIGONOMETRIC 
FUNCTION 
GENERATOR 



25 



TUNABLE OSCILLATOR 
2wn 



CLOCK PULSE GENERATOR 
fc 



/ 



27 



-24 



TUNABLE OSCILLATOR 



L, 



'26 



23 



7 




•19' 



4 



SUBSTITUTE SHEET 



wo 88/02912 



PCT/GB87/00706 




SUBSTITUTE SHEET 



wo 88/02912 



PCT/GB87/00706 



71 



X2(n)= cos2a>Qn 



517 

FIG. 6. 

70 

L 



2 POINT fIR FILTER 



1/^2 .73 



71 



Xj^{n)= cos ff^Qfi 



70 



2 POINT FIR FILTER 




ym'") OUTPUT TO 



T SECONDARY 

SOURCE 



COMPLEX COEFFICIENTS 
OF. EACH HARMONIC AT 
EACH SENSOR 





K 


ADAPTIVE 
ALGORITHM 


i 




f 

TR 


^JURIER 
ANSFORM 









> 


INVERSE 
FOURIER 
TRANSFORM 









OUTPUT WAVEFORMS 
TO EACH SECONDARY 
SOURCE lyn,(n)) 



FIG. 7. 



V/AVEFORMS FROM ERROR SENSORS lei{n)) 

SUBSTITUTE SHEET 



wo 88/02912 



PCT/GB87/00706 





SUBSTITUTE SHEET 



wo 88/02912 



PCT/GB87/00706 



7/7 

FIG. 9. 



100 



10 



5: 




90 ! 



91' 



I 



5 

_^ L 




90 



\ 



\ 





ENGINE 




15; 


1 


16 

\L 


REFERENCE 

SIGNAL 
GENERATOR 


u 





—2 



■13 



ADAPTIVE 
FILTER 



ADAPTIVE 
FILTER 



3 



L 



3 



J 



SUBSTITUTE SHEET 



INTERNATIONAL SEARCH REPORT 

International Apslication No PCT/GB 87 /00706 



I. CLASSIFICATION OF SUBJECT MATTER (if «eveMl classification tympoti apply, indicata all) ' 
AccorcJing to Inrarnational Patent Claaaification (IPC) or to ftotn National Claiaincation and IPC 



IPC 



G 10 K 11/16 



Minimum Documentation Searched ' 



Classification System | 



Classification Symbols 



IPC 



G 10 K 



Documentation Searched other than Minimum Documentation 
to the £it«nt that such Documents are Included In the Fields Searched • 



III. DOCUMENTS CONSIOeRgP TO Bg RSLSVANT* 



Category * | Citation of Document, ii with Indication^ wher» aopfoprlate, oi the relevant passages 




US, A, 4232381 (RENNICK et al . ) 
4 November 1980 

see column 1, line 49 - column 2, 
line 7 

GB, A, 2149614 (SECRETARY OF STATE FOR 

DEFENCE) 12 June 1985 

see page 2, lines 18-62; figure 2 
cited in the application 

GB, A, 782794 (GENERAL ELECTRIC CO.) 
11 September 1957 

see page 2, line 96 - page 3, line 
2; figure 2 

FR, A, 2351466 (SOUND ATTENUATORS LTD) 

9 December 1977 

see claim 1 
cited in the application 
& GB, A, 1577322 



GB, A, 2122052 (THE PLESSEY 
4 January 1984 



CO. ) 



/ 



* Spscial categoriss of citsd docum«nU: t« 
"A" document defining th« ganaraf stat* of th« art which Is not 
cpnsidered to bo oi particular rsUvanca 

ffr^**^ documant but published on or sftsr th« IntsmationaJ 
filing dsts 

"L* document which may throw doubt* on priority cUlmCs) or 
which Is cited to establish tho publication data of anothar 
citation or other epactal raaaen (as spaciflsd) 

"O" document refarring to an oraJ diaclosura. usa, aahlbltlon or 
other means 

"P" documant published prior to tha IntamatlonaJ flUng data but 
later than tha priority data ctalmad 



publlahad after tha International filing data 
?it-H conflict with tha application but 

Invention prtncipla or thaory underiying tha 

iff-U'C* °* Particular ralavanca: tha clalmad Invention 
fl^rH 1 " «Oh«»d«rad nova! or cannot b« considarad to 
Involve an Invantiva step 

ffflf^^r* ®* Dsrtlcular ralavanca;' tha clalmad invention 
cannot be conaidarad to Involva an Invantiva step wnan tha 
oocumant la combinad with one or mora othar euch docu- 
fn'tha'art!^" co'"binatlon baing obvloua to a parson skilled 

"4* documant mambar of tha same patent family 



IV. CIRTIWCATtOW 



Date of tha Actual CompleUon of. tha tntamatlonaJ Saarcft 

8th January 1988 




Data 



of Malting of thla Intamatlonai Saareh Raport 

12 FEB ^388 



41 



/ 



1 -501 34 i 



International Application No. 



PCT/GB 87/0070 



III. DOCUMENTS CONSIDKRCD TO BE KCLIVANT (CONTINUED FROM THE SECOND SHEET) 



Category • ,. Citation of Document, with indicatwn. wT>«fe tporooriate. ol th« r»t«v«nt passages j Ralavant to Claim No 



see abstract 
cited in the application 



Form PCX ISA 210 (aitra shaat) (January 1989) 




ANNEX TO THE INTERNATIONAL SEARCH REPORT 

ON INTERNATIONAL PATENT APPLICATION NO. qB 8700706 

SA 18992 

Tllif m!!mh'l!i",^ the patent family members relating to the patent doeuments cited in the above-mentioned international search reoort 
The members are as contained m the F.uropcan Patent OfTice EDP Hie on '9/01/88 scarcn report. 

The European Patent omce is in no way liable for these particulars nhieh 'arc merely given for the purpose of information. 



Patent document 
cited in search report 


Publication 
date 


Patent family 
member(s) 


Publication 
date 


US-A- 4232381 


04-11-80 


None 






G8-A- 2149614 


12-06-85 


None 






G8-A- 782794 




None 






FR-A- 2351466 


09-12-77 


BE-A- 


854547 


01-09-77 






DE-A,C 


2721754 


24-11-77 






AU-A- 


2483877 


09-11-78 






US-A- 


4153815 


08-05-79 






AU-B- 


507688 


21-02-80 






GB-A- 


1577322 


22-10-80 






SE-A- 


7705504 


14-11-77 






SE-B- 


447937 


22-12-86 


GB-A- 2122052 


04-01t84 


None 







more details about this annex : sec OfTlcial Journal of the Ruropean Patent omce, No. 12/82 



